Magnetic resonance in the cuprates — exciton, plasmon, or tt— mode 
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We re-analyzed the issue whether the resonance peak observed in neutron scattering experiments 
on the cuprates is an exciton, a tt— resonance, or a magnetic plasmon. We considered a toy model 
with on-cite Hubbard U and nearest-neighbor interactions in both charge and spin channels. We 
found that the resonance is predominantly an exciton, even if magnetic interaction is absent and 
d— wave pairing originates from attractive density-density interaction. Our results indicate that one 
cannot distinguish between spin and charge-mediated pairing interactions by just looking at the 
resonance peak in the dynamic spin susceptibility. 



Introduction. The origin of the (tt, tt) spin resonance 
in the cuprates continue to attract interest of the high- 
Tc community. The resonance has been observed in 
four different classes of high- Tc compounds - YBCO, 
Sz2212, T/2201 and i/.gl201-, and the doping variation 
of its energy follows closely the doping dependence of Tc. 
Magnetic resonances have been also recently observed in 
heavy-fermion materials^ and in i^e-pnictides^. 

It is widely accepted that the resonance at Q = (tt, tt) 
is a feedback from the opening of a d— wave pairing gap in 
the fermionic spectrum There is no consensus, however, 
about the driving force. The simplest and most trans- 
parent idea put forward by various groups^ is that the 
neutron resonance is a spin exciton, that is, a resonance 
mode in the spin response function, which emerges due to 
an attractive residual spin interaction between quasipar- 
ticles in a d— wave superconductor. To obtain this mode, 
one can either compute the susceptibility within the low- 
energy model with spin interaction only (no charge com- 
ponent)^, or just calculate the spin susceptibility within 
a conventional RPA for the underlying Hubbard model^ 
- either way one obtains a 5— functional excitonic peak 
at a finite frequency below 2A, where A is the amplitude 
of gap at "hot spots" - kp points separated by Q. 

This simple approach, however, is incomplete as it 
neglects the fact that in a wave superconductor the 
staggered particle-hole, charge spin variable Sq ~ 

{I/2VN) J2k ^I.a<^ap(^k+Q,i3 IS mixcd with the staggered 
d— wave particle-particle charge ±2 variables ttq and 



(tt")*, where tt^ = il/VN)J2kdkCkA'^''cr^)afjCk+Q,fj, 
and dfe = coskx — cosky (Refs. [rll^fioj). Diagram- 
matically, mixed < Sir > response function is given by 
dkGkFk+Q bubbles made out of normal (G) and anoma- 
lous [F) Green's functions^. Such terms are finite in a 
d— wave superconductor at a; ^ 0. 

Because spin and tt responses are coupled, the full spin 
response function is obtained by solving the full 3x3 
set of coupled generalized RPA equations for < SS > , < 
TT, TT > and < Sir > correlators (Refj^). As a consequence, 
the resonance mode emerges simultaneously in spin and tt 



channels, and for the case when both s and tt resonances 
are present, its location lo ~ ujres is in general the solution 
of 



(1) 



where oc (w — ujexc) and x^^ (ui — uj^r) are inverse 
RPA susceptibilities in s and tt channels, each resonating 
at its own frequency, and uj C^^ is the mixing GF terroii. 
If C'i^ = 0, s and tt channels are decoupled, and the spin 
and TT resonances occur at cuexc and , respectively, and 
is are not affected by each other. In general, however, 
the resonance frequency ujres is the solution of ([T]), and 
the full spin and tt susceptibilities near the resonance are 
given by Xs = Zs/{uj - Ures), X-n = Z.„/{uj - ujres) [we 
normalize Z such that for C^j = 0, Zg = l^Z-j^ = at 



: UJe 



and Ztt 



1 , Z„ — at cij 



Eq. (m shows that, in general, there are three possi- 
bilities for the neutron resonance. It can be an exciton, 
which is the case when ujres ~ ^exc and >> Z^r (Refs. 
lii). It can also be tt resonance^, which holds when 
uJres ~ Wtt, and >> Zg. And, finally, it can be a mag- 
netic plasmoE^iS, which is the case when Xs{^) and Xn{^) 
weakly depend on frequency, and the resonance emerges 
due to the mixing between the two channels. In this last 
case, UJres ~ (x7H0)A:ir^(0)/Co)^^^ and is generally dif- 
ferent from both LOexc and o;^. 

Another issue which we consider is the relation be- 
tween the spin resonance and the "glue" for a d— wave 
superconductivity, at least at and above optimal doping, 
where the system falls into moderate coupling regime^^. 
The pairing can be magnetically-mediated^, or it can be 
mediated by a d— wave attraction in the charge channeli^. 
We will analyze whether the location and the residue of 
the resonance can distinguish between the two cases. 

We follow earlier works, use BCS approximation, 
and model the attractive spin-dependent interaction by 
nearest- neighbor Heisenberg exchange term J > 0, 
(Ref.— ), and a d— wave interaction in the charge chan- 
nel by nearest-neighbor density-density interaction V 
(RefiiS). For a repulsive charge interaction {V > 0), 
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there is no tt resonance (i.e., no pole in x-,^ below 2 A) — , 
hence the neutron resonance can be either an exciton or 
a plasmon. For negative 1/, both Xs and Xtt have poles 
below 2 A, and the resonance can be an exciton, a tt- 
resonance, or a plasmon. To distinguish between them, 
we will solve the full 3x3 matrix equation for compare 
the residues in spin and tt channels, (this should show 
whether the resonance is an exciton or a tt resonance), 
and also compare the location of the pole with lo^xc (this 
should show whether or not the resonance is a plasmon) . 
We follow earlier worki^ and require that the value of a 
d— wave gap should agree with ARPES experiments^^ 

Our results show that, to a surprisingly good accuracy, 
the resonance remains an exciton no matter whether the 
pairing is in the spin or in the charge channel. For both 
cases, we found that neither tt— resonance nor the mix- 
ture between s and tt channels affect the location and 
the residue of the resonance residue in any substantial 
way, although the corrections due to mixture are larger 
for the case \y\ >> J. Furthermore, we find that Zg 
can be large enough and the resonance frequency can be 
the experimental AOmeV without placing the system too 
close to an antiferromagnetic instability. 

Our results disagree with the idea about the domi- 
nance of the TT— resonance^, and also somewhat disagree 
with the recent studyi^ which associated the resonance 
with a plasmon rather than an exciton. That work also 
found that the study of the resonance can distinguish be- 
tween spin and charge mechanisms in favor of the former. 
We, on the contrary, found that the resonance is only 
weakly sensitive to the form of the pairing glue. Still, 
we find; in agreement with Ref. [lO|, that the mixing 
between spin and tt channels is not completely negligible 
and has to be taking into account in quantitative studies 
of the cuprates. We caution that, for the charge-mediated 
pairing, the results stronly depend on the magnitude of 
nearest- neighbor attraction V. We have chosen \ V\ which 
yields a BCS gap of i^meV. This \V\ turns out to be 
too small to give rise to tt resonance. For larger \V\, the 
structure of Xs(^) '^iU differ more from an exciton. 

The model. We consider the same model as Lee et al 
(Ref. [Tol ). with on-site Hubbard repulsion and nearest- 
neighbor density-density and spin-spin interactions. 
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FIG. 1: The bare susceptibilities in the spin and tt channels, 
Xs{^) and xii'^) ~ Xiiy respectively (in units of 1/t) The 
susceptibility for the amplitude tt channel X22 is smaller and 
featureless. Note the difference in vertical scales - xSr is larger 
than Xs- 



cosky)/2 is determined from the standard equation 
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9l 



(27r)2 ^(^k - m)2 + A^ff^ 



1. 



(3) 



H y](ek-AOaLaka'^ 



-E 



where ek = — 2t(cos k^ + cos ky) 



(2) 



it' cos fcx COS ky , Hi, 



and Si = [l/2)c\^aiCicr are the particle and spin 
operators on site i (each interaction is counted once). 
A similar model but without U term has been earlier 
considered by Norman and one of us^. 

The soft modes of the system are singlet pairs on 



where V^ = V - 3J/4. Choosing x = 0.12 {^l = -0.94i), 
t' /t = —0.3 and t = 0.433 eV to match the observed 
shape of the Fermi surface and the nodal Fermi veloc- 
ity—, and setting the maximum gap to be A = 357716^, 
we find = V — 3 J/4 = — 0.60t, in agreement with—. 

Generalized RPA susceptibilities We use and A 
as inputs and compute dynamic susceptibilities within a 
generalized RPA scheme which, wc remind, takes into ac- 
cunt the fact that a particle-hole and a particle-particle 
channels are mixed in the presence of a charged conden- 
sate of Cooper pairs. The derivation of the generalized 
RPA equations is rather straightforward and has been 
described hefore^^^. Because of SU{2) spin symme- 
try, it is sufficient to probe only one spin component, 
e.g., restrict, in momentum space with ^0(9) = 
Ai{q) = {TTy{q) + in^{q))/2, and ^2(9) = Al{q). These 
three operators create bosonic excitations with the same 
momentum and spin Sz = I, but with different charges, 
and ±2, respectively. For definiteness, we restrict with 
antiferromagnetic q = Q ~ (7r,7r). 

Generalized RPA equations relate bare and full suscep- 
tibilities: 



(4) 



nearest-neighbor bonds, 



icr^^flj^, spin fluctu- 



ations Sij = {l/2)al^<Taf3aji3, and triplet pairs tt^ 
ata{^cry)af3ajf3. The gap Ak = Ag^ with = (cos A:^; 



where a, 6 = 0,1,2, and Xabi'^) linear response 
functions for the noninteracting system - the Fourier 
transforms of — i0(i)([Aa(t), ^^5(0)]), where the averag- 
ing is over free fermion ground state. Diagrammati- 
cally, Xaai^) the convolutions of GG and FF terms. 
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FIG. 2: The resonance positions for V = and J — 0.8t for 
different U. tufuu is the solution of the full 3x3 set, ^exc IS 
the energy of a spin exciton. The blue dashed line is the edge 
of two particle continuum. We verified analytically that for 
Vi2 > (this case), uifuii > ui^^c- 



while non-diagonal Xabi'^) ^^'"^ terms. All 9 ele- 
ments of Xa bi^) nonzero, but non-diagonal terms 
Xoi and Xi2 vanish at zero frequency, and Xo2 nonzero 
only because of the curvature of the fermionic dispersion 
near kp. It is convenient to rotate the basis in tt, tt* 
plane and introduce, instead of Ai^2, ^1,2(^1 T ^2)/^/^ 
(same notations have been used in earlier works— i^^). In 
this basis, the interaction matrix Fq {, is diagonal due 
to charge conservation, and its diagonal elements are 
Vo ^ -U ~ 2J,Vi , = (y + J/4)/2, i.e., without mix- 
ing Xs = x2/(i - ^0X°), X. - - Vi.2xl)- 

In Fig. [l]we show real and imaginary parts of the bare 
X°(w) = Xoo(^^) and x°(w) = Xn(w) (a phase mode of 
the TT field). The bare X22 smaller than Xii- Note 
that in our notations, staticXs negative. Sharp fea- 
tures at ±0.16t ( ±70meV) are 2A effects, the features 
at higher energies are Van Hove singularities. 

The Results. We present the results for the two ex- 
treme cases cases V = and J = 0. In the first case, 
d— wave superconductivity is magnetically-mediated and 
V^, — —3 J/4. In the second it emerges due to an at- 
traction in the charge channel, and = V. For both 
cases, we used on-site Hubbard U as an extra parame- 
ter that drives the system towards an antiferromagnetic 
instability and brings the resonance frequency down. 

For the first case, Vq = -{U + 1.6t), Vi,2 = O.lt > 0, 
and TT— susceptibility taken alone only develops an anti- 
resonance above the upper edge of two-hole continuum^. 
The issue for this case is whether the resonance is an ex- 
citon or a plasmon. In Fig. [2] we present the results 
of our calculations of the resonance frequency using the 
full 3x3 set and compare them with the RPA result 
for spin-only channel. We see that the energies match 
nearly perfectly, which we believe is a strong indication 
that the resonance is indeed the exciton. On a more care- 
ful look, we found that a near-perfect match is not the 
consequence of the small mixing amplitude {C^ in ([1]) is 
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FIG. 3: The residues Zs, Ztt from the full 3x3 set, and the 
residue of an exciton Zexc for different U. 



roughly 1.3t) but rather the consequence of the fact that 
uj = 50meV is only 0.12i. In Fig. [3]we plot the residues 
of the resonance in spin and tt channels, Zg and Z^^, re- 
spectively, together with the residue of a pure exciton, 
Zexc, which we obtained by eliminating the mixing be- 
tween spin and tt channels. We used a finite broadening 
7 — 0.002 which explains why Z < 1 even for the case of 
a pure exciton. We see that for all U, the residue of the 
resonance is much larger in the spin channel than in the 
TT channel. If the resonance was a plasmon, the residue 
in the spin and tt channels would be comparable. As an 
independent check, we solved a 3 x 3 set with diagonal 
Xaai^) replaced by their static values. The solution in 
this case would be a plasmon (see Eq. ([1])), but we didn't 
find a resonance. 

Fig. [2] also shows that the resonance shifts down from 
2 A (and becomes strong) when |Vo| exceeds roughly 80% 
of the critical |Vo| ~ 2.38t. beyond which antiferromag- 
netic order emerges. Using ^ oc (^/a)^, where ^ is the 
correlation length and a is interatomic spacing, we find 
that this corresponds to ^ ^ 2.5a, while ^ ^ 3.7a is nec- 
essary for the resonance frequency tOres to be AOmeV. 

The results for the second case, J = 0,V = —0.6t, i.e., 
Vb = —U, Vi^2 = — 0.3t, are presented in Figs. |3]and [51 
Now TT-channel becomes are attractive, and exciton, plas- 
mon, and TT resonance are all competing for the dominant 
contribution to the resonance in the full spin susceptibil- 
ity. The three sets of points in Fig. 2] are the solution of 
the full 3x3 set and two approximate sets in which we 
(i) considered the spin channel only (the resonance is an 
exciton) and (ii) approximated diagonal Xaai'^) by their 
static values (the resonance is a plasmon). We see that 
the position of the actual resonance (the full solution) is 
rather close to the position of the spin exciton and the 
two follow the same trend with U, although there is a 
clearly visible difference of about 5 — 10%. The plasmon 
has different dependence on U, and is located below 2A 
only in a narrow range of U. Note that LOfuii < '^exa'^ph 
as it indeed should be (see Eq. ([T|)). In Fig. Owe show 
the residues of the spin and tt components of the full sus- 
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FIG. 4: The energy of the resonance from the solution of the 
fuU 3x3 set (uifuii), and the energies of an exciton and a 
lasmon {uj,;xc and uipi) for the case of nearest-neighbor charge 
interaction J = 0,V = —0.6t as functions of U (see text). The 
blue dashed line is the edge of the two particle continuum. In 
this case, V12 < 0, and ujfuii < ijJres- 
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FIG. 5; The residues Zs, Z^, and Zexc for the case J = 0,V 
-0.6t for different U. 



ceptibilities near iVres and compare them with the residue 
of a pure spin exciton (the case (i) above). We see that, 
when the resonance shifts below 2A and becomes mea- 



surable, its residue in the spin channel is larger than in 
the TT channel and practically coincides with the residue 
of an exciton. Note that the same ^/a 3.7 as in the 
first case is required for the resonance to be at AOmeV. 

These results imply that, even if the pairing is due 
to an attractive nearest neighbor density-density inter- 
action, the resonance in the spin susceptibility still has 
predominantly excitonic character. We caution, however, 
the absence of a substantial tt component of the neu- 
tron resonance is as a numerical rather than a funda- 
mental effect. Namely, for V — ~0.6eV extracted from 
the gap equation, tt— resonance in the absence of broade- 
ing is close to 2A, and a small broadening washes it out. 
We also note that a plasmon does exist in this case, in 
agreement with Ref. — , but in a narrow range of U near 
an antiferromagnetic instability. 

We also performed calculations for several J and V in 
between the two limits, still keeping = V — 3 J/4 = 
— 0.6< to match the gap value. For all cases we found that 
the resonance is predominantly an exciton. The situation 
canges if we abandon BCS gap equation and take larger 
\V\. The larger \V\ is, the stroner the resonance differs 
the exciton. 

To conclude, in this paper we re-analyzed whether the 
resonance peak observed in neutron scattering experi- 
ments on the cuprates is an exciton, a tt— resonance, or 
a magnetic plasmon. We considered a model with on- 
cite Hubbard U and nearest-neighbor interaction in both 
charge and spin channels and found that the resonance is 
predominantly an exciton, even if d— wave pairing orig- 
inates from attractive density-density interaction rather 
than spin-spin interaction. Our results indicate that one 
cannot distinguish between spin and charge- -mediated 
pairings by looking at the resonance peak in the spin sus- 
ceptibility. Other probes like e.g., dispersion anomalies^^ 
or Raman scatteringii are more useful in this regard. 
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